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Abstract. The leptonic tensor for the process e + e~ — > 7 + 7*, which describes the next-to-leading order 
virtual and soft QED corrections to initial state radiation in e + e~ annihilation with emission of an extra 
virtual photon decaying into hadrons, is calculated. A Monte Carlo generator for the reaction e + e~ — > 
7 + 7r + 7r~ has been set up which includes these corrections. It thus describes configurations where the 
invariant mass of the hadrons plus photon is very close to y/s. Predictions for cms energies of 1 to 10 GeV, 
corresponding to the energies of DAPHNE and B-meson factories, are presented. The possibility for an 
accurate measurement using tagged photons of a(e + e~ —> hadrons), which plays an important role in the 
theoretical description of the muon anomalous magnetic moment and the running of the electromagnetic 
coupling, is discussed. 
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1 Introduction 

Electroweak precision measurements have become one of 
the central issues in present particle physics. The recent 
measurement of the muon anomalous magnetic moment 
(g — 2) M reported in shows a discrepancy at the 2.6 
standard deviation level with respect to the theoretical 
evaluation of the same quantity [Q. This disagreement, 
which has been taken as an indication of new physics, has 
triggered a raving and somehow controversial deluge of 
non-Standard Model explanations. 

One of the main ingredients of the theoretical pre- 
diction for the muon anomalous magnetic moment is the 
hadronic vacuum polarization contribution || which more- 
over is responsible for the bulk of the theoretical error. It 
is in turn related via dispersion relations to the cross sec- 
tion for electron-positron annihilation into hadrons, o~had- 
This quantity plays also an important role in the evolu- 
tion of the electromagnetic coupling chqed from the low 
energy Thompson limit to high energies |,|. This indeed 
means that the interpretation of improved measurements 
at high energy colliders like LEP, LHC, or Tevatron de- 
pends significantly on the precise knowledge of (Jhad- 

The evaluation of the hadronic vacuum polarization 
contribution to the muon anomalous magnetic moment, 
and even more so to the running of cxqed, requires the 
measurement of (Jhad over a wide range of energies. Of par- 
ticular importance is the low energy region around 2 GeV, 
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where o~had is at present experimentally poorly determined 
and only marginally consistent with the predictions based 
on pQCD. New efforts are therefore mandatory in this 
direction, which could help to either remove or sharpen 
the discrepancy between theoretical prediction and exper- 
imental results for (g — 2) M . 

The feasibility of using tagged photons at high lumi- 
nosity electron-positron storage rings, like the 0-factory 
DAPHNE or at -B-factories, to measure a tad has been 
proposed and studied in detail in |||] (See also In 
this case, the machine is operating at a fixed cms energy 
and initial state radiation is used to reduce the effective 
energy and thus the invariant mass of the hadronic sys- 
tem. The measurement of the tagged photon energy helps 
to constrain the kinematics. Preliminary experimental re- 
sults using this method have been presented recently by 
the KLOE collaboration M. Large event rates were also 
observed by the BaBar collaboration . 

In this paper we consider the next-to-leading order 
(NLO) QED corrections to initial state radiation (ISR) 
in the annihilation process e + e~ — > 7 + hadrons. After 
factorizing the hadronic contribution, the leptonic tensor, 
which contains virtual and soft photon corrections, is cal- 
culated. An improved Monte Carlo generator including 
these terms has been set-up. Predictions for the exclusive 



channel e + e 



"7 at cms energies of 1 to 10 GeV, 



corresponding to the energies of DAPHNE and £?-meson 
factories, are presented. The comparison with the EVA |j| 
Monte Carlo, simulating the same process at leading-order 
(LO), is performed. 
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Fig. 1. Initial state radiation in the annihilation process 
e + e~ — > 7 + hadrons at the Born level. 
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It is symmetric under the exchange of the electron and the 
positron momenta. Expressing the bilinear products yt by 
the photon emission angle in the center of mass frame 



2/1,2 



l-q 2 



(1 =F (3 cos < 



/3 = \/l - 4m 2 , 



and rewriting the two-body phase space 



2 The leptonic tensor at leading order 

Consider the e + e~ annihilation process 

e+(p 1 ) + e-(p 2 )^ 1 *(Q)+ 1 (k 1 ) , 

I 

hadrons 



(1) 



where the virtual photon decays into a hadronic final state 
and the real one is emitted from the initial state (figure |I|) . 
The differential rate can be cast into the product of a lep- 
tonic and a hadronic tensor and the corresponding factor- 
ized phase space 



da = Y L »» H ' 1 ' / d&2(pi,P2;Q,ki)d$n(Q;qi,;<ln)-^r 



dQ 2 
2tt 



(2) 

where d<P n (Q', Qi,',Qn) denotes the hadronic n-body phase- 
space with all the statistical factors coming from the hadro- 
nic final state included. 

For an arbitrary hadronic final state, the matrix ele- 
ment for the diagrams in figure ^ is given by 



M = M % - 



1 VWHfi-yi+meh" 



(3) 
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where is the hadronic current. Summing over the po- 
larizations of the final real photon, averaging over the po- 
larizations of the initial e + e~ state, and using current con- 
servation, Q M J p = 0, the leptonic tensor 



jliv 



can be written in the following form: 

(Ana/s) 2 \ (2m 2 q 2 (\- q 2 ) 2 2q 2 + y\ + y\ 
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it is evident that expression (|]) contains several singular- 
ities: soft singularities for q 2 — > 1 and collinear singular- 
ities for cos# — > ±1. The former are avoided by requir- 
ing a minimal photon energy. The later are regulated by 
the electron mass. For s m 2 , the expression (0) can 
be nevertheless safely taken in the limit m e — > if the 
emitted real photon lies far from the collinear region. In 
general, however, one encounters spurious singularities in 
the phase space integrations if powers of m 2 = m 2 /s are 
prematurely neglected. 

The physics of the hadronic system, whose description 
is model dependent, enters only through the hadronic ten- 



(7) 



where the hadronic current has to be parametrized through 
form factors. For two charged pions in the final state, the 
current 



(8) 



where q^+ and q w - are the momenta of the 7r + and tt~ 
respectively, is determined by only one function, the pion 
form factor F 2 ^ pd[ |. The hadronic current for four pions 
exhibits a more complicated structure and has been dis- 
cussed in p^ |. 

After integrating the hadronic tensor over the hadronic 
phase space, one gets 

J H^ n (Q; qi ,;q n ) = ^{Q»Q V - g^Q 2 )R{Q 2 ) , 

(9) 

where R{Q 2 ) = a(e + e~ — > hadrons) / o(e + e~ — > 

After the additional integration over the photon angles, 

the differential distribution 



dQ z Ss [s(s — Q z ) 



(10) 



is obtained. If instead the photon polar angle is restricted 
to be in the range 9 m in < < n — 6. 



this differential 
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Fig. 2. One- loop corrections to initial state radiation in 
e + e~ — > 7 + hadrons. 



distribution is given by 
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In the later case, the electron mass can be taken equal 
to zero before integration, since the collinear region is ex- 
cluded by the angular cut. The contribution of the two 
pion exclusive channel can be calculated from eq.([lOj) and 
eq.(jll]) by the substitution R(Q 2 ) -> (1 - 4m 2 /Q 2 ) 3 / 2 
I F 2 „(Q 2 ) | 2 /4. 



3 Virtual and soft photon corrections to the 
leptonic tensor 

At NLO, the leptonic tensor receives contributions both 
from one- loop corrections (figure |^) arising from the inser- 
tion of virtual photon lines in the tree diagrams of figure [l] 
and from the emission of an extra real photon from the 
initial state. In this paper, we consider only the emission 
of soft photons. The contribution of a second hard photon 
will be considered in a separate work [jl3| . 

The one-loop matrix elements contribute to the lep- 
tonic tensor through their interference with the lowest or- 
der diagrams (figure [l]) . They contain ultraviolet (UV) and 
infrared (IR) divergences which are regularized using Di- 
mensional Regularization in D = 4 — 2e dimensions. The 
UV divergences are renormalized in the on-shell scheme. 
The IR divergences are canceled by adding the contribu- 
tion of an extra soft photon emitted from the initial state 
and integrated in the phase space up to an energy cut- 
off Ery < Wy/s far below y/s. The result, which is finite, 
depends, however, on this soft photon cutoff. Only the 
contribution from hard photons with energy E 1 > w^/s 
would cancel this dependence. 

The algebraic manipulations have been carried out with 
the help of the FeynCalc Mathematica package Us- 
ing standard techniques |l5[ ] , it automatically reduces the 
evaluation of the one-loop contribution to the calculation 
of a few scalar one-loop integrals, listed in appendix |b], 
and performs the Dirac algebra. 



At NLO, the leptonic tensor has the following general 
form: 
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The scalar coefficients et^ and a_i allow the following ex- 
pansion 



(0) . OL (1) 



0-1 = - a a ] . (13) 

7T 



The LO coefficients a-^ can be directly read from eq.([|). 
For vanishing electron mass 



vl) 



= 0, 



,(()) 
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(14) 



Below, we use m 2 < s,Q 2 and neglect terms proportional 
to m 2 , which is a valid approximation if the observed pho- 
ton is far from the collinear region. The imaginary anti- 
symmetric piece proportional to a_i appears for the first 
time at the NLO. The leptonic tensor remains therefore 
fully symmetric only at the LO. 

After adding the one-loop and the soft contribution, 
we end up with expressions for the NLO coefficients a^' 

and a_\ which can be found in appendix [A|. As a test of 
our calculation, the leptonic tensor has been contracted 
with (q^q u — q 2 g tiu ). The result: 



L^iq^qv - q 2 g^) = (q^q v - q 2 g liv ) 



(15) 



-log(4u; 2 )[l + log(m 2 )] 

2j/i?/2 



1 - 



ym 



1 



22/12/2 

y\ - 2y 2 _ 

2(1 -y 2 ) 4(1 -y 2 ) 2 
1 - 



2<z 2 + y\ + y\ 



log(^) 



1 



2/2 



2(1 -g 2 ) 4(1 - Vl ) 



where 



LQ V {q^qu - q 2 9nv) = {Ana/sf 



2{2q 2 



(yi 



V2. 



vl) 



q 2/12/2 



(16) 



reproduces, up to coupling constants and global factors, 
the one given in |l6| for the QED radiative corrections to 
the reaction e + e~ — > Z7 . 

The leptonic tensor can be cast, from eq.([l2]) and ap- 
pendix into the form: 



1 



• log(m 2 ) - 2 



log(4w 2 )[l + log(m 2 )] 
.2 - 
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Fig. 3. Initial state radiation in the process e + e~ — > 7r + 7r~7 
at the Born (blue) and NLO level for yfs = 1.02 GeV. Virtual 
plus soft contributions are considered with a soft cutoff equal 
to w — 0.01 at fixed order (yellow) and exponentiated (red). 
The final state radiation contribution (including interference 
with ISR) is also shown (green). The cuts are: 7° < 6* 7 < 20° 
or 160° < By < 173°, 30° < < 150° and the energy of the 
observed photon _E 7 > 0.02 GeV. 
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Fig. 4. Initial state radiation in the process e + e — > tt + -k 7 
at the Born (blue) and NLO level for y/s = 10.6 GeV. Virtual 
plus soft contributions are considered with a soft cutoff equal 
to w = 0.001 at fixed order (yellow) and exponentiated (red). 
The cuts are: 7° < 6 y < 173° and the energy of the observed 
photon £ 7 > 0.02 GeV. 



where the expected soft and collinear behaviour is man- 
ifest. The first term, log(4u; 2 )[l + log(m 2 )], coming from 
the emission of soft photons, gives a large contribution 
which eventually could spoil the improvement expected 
from a NLO analysis. The soft cutoff w should be small 
enough to justify the soft photon approximation. On the 
other hand, a very small value of w would result in a mean- 
ingless negative cross-section. To overcome this difficulty, 
the contribution of a second hard photon, with energy 
E 7 > w-y/s, should be added [|l3], thus canceling the in- 
dependence. However, one may also limit the analysis to 
configurations where the additional radiated photon is re- 
quired to be soft, by constraining the invariant mass of 
the hadron + photon system. For small w the following 
exponentiation can be used: 



1 - - log(4u; 2 )[l + log(m 2 )] -> (2iu) 

7T 



^[l+log(m 2 )] 



(18) 



This accounts for the leading soft logarithms to all orders 
in perturbation theory fTij| and vanishes in the limit w — > 
0, as expected. 

Another contribution which can be rather large and 
was not considered up to now, is the insertion of vacuum 
polarization corrections to the virtual photon line in fig- 
ure 0. Its inclusion does not affect the other features of 
our calculation. It introduces a correction which is pro- 
portional to the Born leptonic tensor and can be reab- 
sorbed in the choice of a running coupling constant. In 
order to take into account also the higher orders, a factor 
1/(1 — Syp) where Syp is the vacuum polarization con- 
tribution can be included. In the present version of this 
program, this factor has been dropped for simplicity. 

Motivated by these considerations, the following im- 
proved leptonic tensor can be defined: 



(2wY 



[l+log(m 2 )] 



1VP 



(19) 



log(m 2 ) 
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4 Monte Carlo simulation 

A Monte Carlo generator has been built. It simulates the 
process e + e~ — > 7r + 7r~7 where the photon is emitted from 
the initial state, with a large angle with respect to the 
beam direction. It is based on EVA || and includes the 
NLO corrections described in the previous sections [ . The 
program is built in a modular way such that the simulation 
of other exclusive hadronic channels can be easily included 
while keeping the factorization of the initial state QED 
corrections. Configurations with an additional hard pho- 
ton are not yet included. Hence the generator can only be 



1 The default version of the program is based on eqs( pc| - |24] ) 
and is thus valid in the limit m^/s < 1. As an alternative, it is 
also possible to run the program with formulae which include 
the complete m e dependence. 
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used to describe configurations where the invariant mass of 
the hadronic system plus the tagged photon is close to the 
total energy of the collision, M(j+hadrons) > (l — w)y/s. 
A generator which includes also emission of two hard pho- 
tons will be presented in a subsequent publication. 

In figure [5[ the differential distribution in the invari- 
ant mass of the hadronic system is shown for the process 
e+e- -> 7r+7r~7 at DAPHNE energies, y/s = 1.02 GeV. 
In principle, initial and final state radiation (FSR) would 
be required for the complete simulation of the process. 
However, the cuts: pions in the central region, 30° < 8^ < 
150°, and photons close to the beam and well separated 
from the pions, 7° < 7 < 20° or 160° < 7 < 173°, select 
configurations dominated by initial state radiation |^,[|. 
For comparison, the contribution of the final state radi- 
ation at Born level is also shown. The corrections to the 
Born approximation reach up to 20% for a soft cutoff of 
w = 0.01 if the fixed order expression (0) is used, or 17% 
if the exponentiation of the soft cutoff, eq. ( |l8| ) , is applied. 

At B-factories, y/s — 10.6 GeV, the situation is slightly 
different. Because of the higher center of mass energy, 
the hadronic system and the real photon are mainly pro- 
duced back to back. The FSR is therefore already sup- 
pressed without imposing additional kinematical cuts on 
the hadronic system. For these high energies one might 
in addition include photons emitted also at larger an- 
gles. From the analysis of EVA || one finds that the 
events are always dominated by ISR - a consequence of 
the strong form factor suppression of FSR at high ener- 
gies and large angles between pions and photons. Figure [| 
shows the differential distribution in the invariant mass of 
the hadronic system at yfs = 10.6 GeV for photons in the 
region 7° < 7 < 173° and no constraints to the hadronic 
system. The FSR amounts to less than 10~ 3 and is not 
shown. The NLO corrections are again dominated by the 
value of the soft photon cutoff. 

5 Leading log resummation 

From the results of the previous section, it is clear that 
the fixed order correction is dominated by the value of the 
soft photon cutoff w. Even if the exponentiation is applied, 
which accounts for the leading soft logarithms to all or- 
ders in perturbation theory, the situation does not improve 
much. The addition of the contribution of a second hard 
photon would cancel the strong w-dependence, however, 
large logarithms of collinear origin, L — log(s/mg), would 
remain. Because these large logarithms will show up in 
all orders of perturbation theory, resummation techniques 
have been developed which are constructed to include the 
dominant higher order terms. The method of structure 
functions Jl^,[T^] accounts well for these corrections. It 
takes into account all the leading logarithmic corrections 
G(aL), coming from virtual, soft, and hard photon con- 
tributions to all orders in perturbation theory. 

In H these effects were considered in the simulation 
of the process e + e~ — > 7r + 7r~7. The formulas of [^0| were 
used which include terms up to order a 2 L 2 . A minimal in- 
variant mass of the 7r + 7r _ 7 system was required in order 




Q 2 (GeV 2 ) 



Fig. 5. Subleading NLO contribution to initial state radiation 
in the process e + e~ — > 7r + 7r _ 7. Same cuts as in figure ^| 

to reduce the kinematic distorsion of the events due to the 
collinear initial state radiation. For a sufficiently tight cut 
this corresponds to the situation discussed in the present 
paper. For a minimal invariant mass of the 7r + 7r~7 sys- 
tem of 0.9 GeV 2 , a negative 6% correction to the Born 
approximation was found. 

Subleading effects, which correspond to the terms 
of our leptonic tensor in eq.(|l7j), are not taken into account 
within this approximation. In figure |^, the contribution of 
these subleading terms to the differential distribution of 
figure H is shown. They amount to 0.6% of the Born ap- 
proximation close to the p-resonance, where the rate is 
maximal, and to 1% for small values of Q 2 , where, how- 
ever, the cross section rate is small. Similar values are 
found for higher center of mass energies. The matching 
of the resummed result, through the structure function 
technique, with the fixed order result is under study [fl3|| . 

6 Conclusions 

ISR at high luminosity e + e~ colliders (4> and B-factories) 
is an alternative to the direct measurement of o(e+ eT — » 
hadrons) giving access to a wide range of energies, from 
threshold to the center of mass energy of the collider. 
The NLO QED correction to ISR, in the form of a lep- 
tonic tensor, has been calculated and included in a Monte 
Carlo generator, which was compared with EVA || for 
the tt + tt~ exclusive channel. The modular structure of 
the calculation, independent of the final state hadronic 
system, is such that other hadronic channels or improved 
parametrizations of the hadronic current can be easily in- 
cluded. 

From our results it is clear that the fixed order correc- 
tion is dominated by the value of the soft photon cutoff 
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w. Only the contribution of a second hard photon |13| 
would cancel this strong ^-dependence. But even if this 
contribution is added, large logarithms of collinear origin, 
L = log(s/rrig), would still remain. Furthermore, they will 
show up in all orders of perturbation theory. A consis- 
tent way to resum such leading log terms to get her with 
subleading effects is also under consideration [[L3| . 
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A Tensor coefficients 



In this appendix, we collect our results for the scalar coef- 
ficients of the leptonic tensor ([l2]) at the NLO arising from 
virtual and soft photon contributions. Our expressions are 
valid in the limit of small electron mass. 

For the coefficient proportional to g^, we get 
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with 
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q A q z 
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where Li 2 is the Spence or dilogarithm function and y^, 
q 2 and m 2 have been defined in eq.(||). The quantity w 
denotes the dimensionless value of the soft photon energy 
cutoff: E 1 < Wyfs. The coefficient in front of the tensor 



structure PiPi, is given by 
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that in front of p^P 2 can De obtained by symmetric con- 
siderations, exchanging the positron with the electron mo- 
menta 
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For the symmetric tensor structure (PiP 2 + P2P1)' we 
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Finally, the antisymmetric coefficient a_i, accompanying 
(P\P 2 -P 2 P\)i reads 



log(m 2 y 2 ) , 1 - o 2 
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B Loop integrals and soft photon contribution 

The Passarino-Veltman procedure [|l5| allows to reduce 
the calculation of any one-loop amplitude to the evalua- 
tion of a few scalar one-loop integrals. In this appendix, 
we collect the scalar one-loop integrals needed for our cal- 
culation and give their expression in the limit of small 
electron mass. Ultraviolet (UV) and infrared (IR) diver- 
gences appear in the one-loop calculation. Dimensional 
regularization in D = 4 — 2e dimensions is used to regu- 
larize both kind of divergences. The UV divergences are 
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renormalized using the on-shell renormalization scheme 
The soft photon contribution to the leptonic tensor can- 
cels the remaining IR divergences. 

A few two-, three-, and four-point scalar one-loop in- 
tegrals enter our calculation. Expression for the two-point reads 
scalar integrals are simple and well known. The general 
three-point scalar one-loop integral is defined by 



B.l Initial state soft photon emission 



The contribution to the leptonic tensor of a soft photon 
with momentum k and energy E < Wy/~s, with w <C 1, 



rt I 2 / \2 2 2 2 2\ 

Co(PaAPa-Pb) ,p b ,m 1 ,m 2 ,m 3 ) 
d D k 1 



■ 1C 2 4-ZJ 
-ilDTT \l 



^soft - 



-4e 



fi 



i-D 



(27T)D [fc2 _ m 2] p _ pa) 2 _ m 2] [{k _ pb) 2 _ m 2 



2(2tt) d - 1 
Pi P2 



w^fs 



2pi ■ k 2p2 ■ k 



E D - 3 dEdf2 



(30) 



(25) 



Integrated in the soft photon phase space jn| and in the 



Four different three-point scalar one-loop integrals are needed limit of sma11 electron mass > we S et 

L soft = %{i A - 10 S( 4W2 )] t 1 + l0 S( m2 )] 

log 2 (rn 2 ) 



C01 = C (fe-A: 1 ) 2 ,0,77i 2 ,0,TO 2 ,TO 2 ) , 
C02 = Co(m 2 , s, to 2 , 0, to 2 , m 2 ) , 
C03 = C ((Pi - h) 2 , Q 2 ,ml 0, to 2 , to 2 ) 
C04 = C* (Q 2 ,s,0,to 2 ,to 2 ,to 2 ) , 

i = 1,2, and one scalar box 



log(m ) 



(26) 



L 



(31) 



DO = -HQtt 2 ^- d 



d D k 
(2ir) D 



After UV renormalization, the renormalized one-loop ma- 
trix elements Mi of the diagrams in figure |^, which con- 
tribute at the NLO to the leptonic tensor through their 
(27) interference with the Born diagrams, give the following IR 
contribution 



k 2 [(k +pi) 2 - m 2 ][(k + pi - ki) 2 - m 2 ][(k - pj) 2 - m 2 ] 
with j ^ i. 

In the limit of vanishing electron mass: to 2 <C s, Q 2 , Hi, 
the following results are obtained: 



i=l,4 

a 

7T 



Zi [log(m 



1] L; 







(32) 



coi 

C02 = s" 1 
C03 



1 2 TO/ 

2 J/i 



7T" 

Y 

2/ 



fog(TO 2 )Z\ 



log (to 
2"" 

■ Vi 



27T 2 

3 



in A 



bumming up eq. @ and eq.(|32|), a finite result is found. 
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They can be compared, for instance, with the one-loop 
integrals quoted in where a fake photon mass A was 
used to regularize the IR divergences. The identification 
A <-* log(A 2 /s) allows to pass from one scheme to the 
other. 



H. N. Brown et al. [ Muon g-2 Collabo ration], Phys. Rev 



Lett. 86 (2001) 2227 [tiep-ex/0102017 



V. W. Hughes and T. Kinoshita, Rev. Mod. Phys. 71 
(1999) S133. 

W. J. Marciano a nd B. L. Rober ts, friep-ph/0105056 |. 
F. Jegerlehner, |hep-ph/0104304| . M. Davier an d 
A. Hocker, Phys. Lett. B 435 (1998) 427 fcep-ph/9805470 |. 
S. Ei delman and F. Jegerlehner, Z. Phys. C 67 (1995) 
585 [frep-ph/9502298fl . K. Adel and F. J. Yndurain, 
[ hep-ph/9509378^ 
4. H. Burkhardt and B. Pietrzyk, LAPP-EXP-2001-03. 
A. D. Martin, J. Outhwaite and M. G. Ryskin, 



[ liop-ph/001223l| 1; Phys. Lett. B 492 (2000) 69 flhep 



ph/0008076]. J. Erler, Phys. Rev. D 59 (1999) 054008 pep 



ph/980345;]. J. H. Kiihn and M. St einhauser, Phys. Lett 



B 437 (1998) 425 hep-ph/9802241 



5. S. Binner, J H. Kiihn and K . Melnikov, Phys. Lett. B459 
(1999) 279 ||hep-ph/990239S |. 

6. K. Melnikov, F. Nguyen, B. Valeriani and G. Venanzoni, 
Phys. Lett. B477 (2000) 114 jhep-ph/0001064 

7. S. Spagnolo, Eur. Phys. J. C 6 (1999) 637. 

8. V. A. Khoze, M. I. Konchatnij, N. P. Merenkov, 
G. Pancheri, L. Trentadu e and O. N. She khovzova, Eur. 
Phys. J. C 18 (2001) 481 ||hep-ph/0003313 l. 



8 G. Rodrigo, A. Gehrmann, M. Guilleaume and J.H. Kiihn: NLO QED corrections to ISR in e + e annihilation 



9. A. Denig, talk given at Physics at Intermediate Energies 
Workshop, SLA C, April 2001. M . Adinolfi et al. [KLOE 
Collaboration], ftiep-ex/0006036 |. 

10. Solodov, talk given at Physics at Intermediate Energies 
Workshop, SLAC, April 2001. 

11. J. H. Kiihn and A. Santamaria, Z. Phys. C 48 (1990) 445. 

12. H. Czyz and J. H . Kuhn, Eur. Phys. J. C 18 (2001) 497 
hep-ph/0008262 |. 

13. H. Czyz, J. H. Kiihn and G. Rodrigo, in preparation. 

14. R. Mertig, M. Bohm and A. Denner, Comput. Phys. Com- 
mun. 64 (1991) 345. 

15. G. Passarino and M. Veltman, Nucl. Phys. B 160 (1979) 
151. 

16. F. A. Berends, G. J. Burgers and W. L. van Neerven, Phys. 
Lett. B177 (1986) 191. 

17. D. R. Yennie, S. C. Frautschi and H. Suura, Annals Phys. 
13 (1961) 379. F. Bloch and A. Nordsieck, Phys. Rev. 52 
(1937) 54. 

18. E. A. Kuraev and V. S. Fadin, Sov. J. Nucl. Phys. 41 
(1985) 466. G. Altarelli and G. Martinelli, in Physics at 
LEP, CERN- Yellow Report 86-06, eds. J. Ellis and R. 
Peccei (CERN, Geneva, February 1986). O. Nicrosini and 
L. Trentadue, Phys. Lett. B196 (1987) 551. 

19. F. A. Berends, W. L. van Neerven and G. J. Burgers, Nucl. 
Phys. B297 (1988) 429; Erratum Nucl. Phys. B304 (1988) 
921. 

20. M. Caffo, H Czyz and E. R emiddi, Nuovo Cim. A110 
(1997) 515 |iep-ph/9704443C ; Phys. Lett. B327 (1994) 
369. 

21. G. Rodrigo, A. Santamaria an d M. Bilenky, J. Phys. G25 
(1999) 1593 friep-ph/9703"36C ]. 



